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ABSTRACT 


We present the first homogeneous catalog of Kepler, K2, and T'ESS host stars and the corresponding 
catalog of exoplanet properties, which contain 7993 stars and 9324 planets, respectively. We used 
isochrone fitting and Gaia DR3 photometry, parallaxes, and spectrophotometric metallicities to com- 
pute precise, homogeneous Teg, log g, masses, radii, mean stellar densities, luminosities, ages, distances, 
and V-band extinctions for 3248, 565, and 4180 Kepler, K2, and TESS stars, respectively. We compared 
our stellar properties to studies using fundamental and precise constraints, such as interferometry and 
asteroseismology, and find residual scatters of 2.896, 5.696, 5.096, and 3196, with offsets of 0.296, 1.096, 
1.296, and 0.7% between our Teg, radii, masses, and ages and those in the literature, respectively. In 
addition, we compute planet radii, semimajor axes, and incident fluxes for 4281, 676, and 4367 Kepler, 
K2, and TESS planets, respectively, and find that the exoplanet radius gap is less prominent in the 
K2, TESS, and combined samples than it is in the Kepler sample alone. We suspect this difference 
is largely due to heterogeneous planet-to-star radius ratios, shorter time baselines of K2 and T'ESS, 
and smaller sample sizes. Finally, we identify a clear radius inflation trend in our large sample of hot 
Jupiters and find 150 hot sub-Neptunian desert planets, in addition to a population of over 400 young 


host stars as potential opportunities for testing theories of planet formation and evolution. 


1. INTRODUCTION 


'The Kepler, K2, and TESS Missions have enabled the 
detection of thousands of transiting exoplanets, from 
close-in hot Jupiters to Earth-analogs, transforming our 
understanding of exoplanet demographics and planet 
formation and evolution. Until now, most analyses have 
focused either on the detection and characterization of 
individual systems (e.g. Holman et al. 2010; Howell et al. 
2012; Hébrard et al. 2013; Ciceri et al. 2015; Vanderburg 
et al. 2015; Schlieder et al. 2016; Cloutier et al. 2020) or 
the characterization of populations of exoplanets from a 
particular mission (Dressing & Charbonneau 2015; Pe- 
tigura et al. 2017; Johnson et al. 2017; Fulton et al. 
2017; Fulton & Petigura 2018; Berger et al. 2018b; Stas- 
sun et al. 2019; Berger et al. 2020a,b; Hardegree-Ullman 
et al. 2020; Zink et al. 2021), but none have performed 
a homogeneous characterization of the entire population 
of host stars to transiting exoplanets. 
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Homogeneous determinations of stellar properties are 
important because of the potential for strong systemat- 
ics between different observables methodologies convert- 
ing those observables into fundamental physical param- 
eters. For instance, it has been shown that the same 
spectrum analyzed by different pipelines produces dif- 
ferent results (Torres et al. 2012), with discrepancies ex- 
ceeding uncertainties in some cases (Tayar et al. 2022). 
While it is difficult to determine fundamental parame- 
ters accurately in an absolute sense, comparing stellar 
parameters homogeneously and then comparing them to 
each other remains an easier task. 

Here, we present the first homogeneous characteriza- 
tion of the population of Kepler, K2, and TESS transit- 
ing planet host stars using isoclassify (Huber et al. 
2017; Berger et al. 2020a), only leveraging observables 
from Gaia DR3 (Gaia Collaboration et al. 2016, 2021, 
2022; Babusiaux et al. 2022): parallaxes, Ggp and Grp 
photometry, spectrophotometric metallicities, and mea- 
sured positions. We re-derive stellar Tog, logg, radii, 
masses, densities, luminosities, and ages for 7993 Ke- 
pler, K2, and TESS hosts and compare our performance 
to current best estimates from interferometry, astero- 
seismology, and clusters. Finally, we re-compute planet 
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Table 1. CKS-Gaia Spectrophotometric Metallicites 


Bounds Adopted Metallicity 


[Fe/H]adopt = —5.30 + 0.591 [M/H]a 
— 0.0257 [M/H] + 0.00193 Teff,G 
— 0.000000171 T2. 

[Fe/H] adopt = 0.16 + 0.605 [M/H]a 


1. 4000 < Tore,g < 7000 K 
N [M/H]a > -2.0 dex 


2. Other 


NoTE—CKS-based polynomial relations used to fix Gaia spectrophoto- 
metric metallicities. The G denotes Gaia spectrophotometric param- 
eters, while [Fe/H]aaopt represent our adopted metallicities. We plot 
the results of these relations in Figure 1. 


radii, semi-major axes, and incident fluxes by combin- 
ing our new stellar properties with results from planet 
transit-fitting pipelines and compare the host star and 
planet populations of Kepler, K2, and TESS. 


2. METHODOLOGY 
2.1. Sample Selection 


For Kepler hosts, we used the Cumulative KOI table 
at the NASA Exoplanet Archive (Batalha et al. 2013; 
Burke et al. 2014; Rowe et al. 2015; Mullally et al. 2015; 
Coughlin et al. 2016; Thompson et al. 2018; NASA Ex- 
oplanet Archive 2021)!, for K2 hosts, we utilized the 
Zink et al. (2021) catalog, and for T'ESS hosts, we used 
the TESS Project Candidates table (e.g. Guerrero et al. 
2021) from the NASA Exoplanet Archive”. We followed 
a procedure identical to Berger et al. (2018b) to cross- 
match Kepler hosts to Gaia DR3 sources, but instead 
used TOPCAT (Taylor 2005) to perform the sky cross- 
match. For K2and TESS hosts, we first selected all Gaia 
DR3 sources within 4" and determined the minimum of 
the resulting angular distance histograms (1.5" for K2 
and 1.0" for TESS). We then removed all matches with 
angular distances larger than that angular distance. For 
the remaining matches where one Ecliptic Plane Input 
Catalog (EPIC, Huber et al. 2016a) or TESS Input Cat- 
alog (TIC, Stassun et al. 2019) star has multiple Gaia 
DR3 sources, we choose the smallest absolute G — Kp 
and G—T discrepancy relative to the median magnitude 
difference for a star with that Kp or T magnitude flag. 

From there, we removed stars without positive Gaia 
parallaxes or Gpp/Grp photometry. We next removed 
all planet candidates already flagged as false positives 
in each planet sample. After this step, we retained 3309 
Kepler, 577 K2, and 4198 TESS host stars. 


2.2. Input Parameters 
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Figure 1. California-Kepler Survey (CKS, Petigura et al. 
2017) spectroscopic metallicities versus those from Gaia DR3 
spectrophotometry (Gaia Collaboration et al. 2022; Andrae 
et al. 2022; Creevey et al. 2022; Fouesneau et al. 2022). In 
the top panel, we plot a direct comparison between the CKS 
and Gaia metallicities, where the solid red line represents 
agreement and the colored curves show the best-fit polyno- 
mial relations used to correct Gaia metallicities for stars of 
different Ter. The bottom panel shows the residuals between 
the Gaia and CKS data after applying the corrections to the 
Gaia metallicities and the black-dashed line represents agree- 
ment. Row 1 of Table 1 contains the polynomial relation used 
to correct Gaia metallicities. 


'To create a homogeneous catalog for the combined Ke- 
pler+K2+ TESS host star samples, we used Gaia DR3 
Gpp and Grp as our photometric constraints. We chose 
these two passbands due to their mmag precision (Riello 
et al. 2021), color-T.g sensitivity, and Grp’s relative in- 
sensitivity to extinction. The typical host star in our 
sample experiences Ay — 0.15 mag, which we correct 
for with an extinction map. In addition, Gaia pho- 
tometry has a S1" angular resolution (de Bruijne et al. 
2015) compared to 2MASS’s 4” (Skrutskie et al. 2006), 
which minimizes issues of photometric blending. "We 
used Riello et al. (2021) equations C.1-C.3 to correct 
saturated photometry, where relevant. While we did 
not modify the formal, input uncertainties on the Gaia 
photometry, their mmag precision makes them difficult 
to use with a relatively sparse pre-computed grid like 
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isoclassify uses; hence, we will detail the required 
modifications to isoclassify in §2.3. 

In addition to the Ggp and Grp photometry, we uti- 
lized Gaia DR3 parallaxes (Lindegren et al. 2021a) with 
zeropoint corrections (Lindegren et al. 2021b) and J2000 
positions. Where possible, we used spectrophotometric 
metallicities (mh_gspphot, Creevey et al. 2022; Foues- 
neau et al. 2022; Andrae et al. 2022). Given that the 
spectrophotometric metallicities exhibit strong system- 
atic errors (Andrae et al. 2022), we corrected them 
using California-Kepler Survey (CKS, Petigura et al. 
2017) metallicities by fitting a quadratic polynomial as a 
function of spectrophotometric Teg (teff gspphot) and 
metallicity to the CKS metallicities. We only applied 
this relation to stars with Gaia spectrophotometric Tog 
between 4000 and 7000 K and metallicities >—2.0 dex. 
For stars outside that range, we used a line representing 
the median offset between the Gaia spectrophotomet- 
ric metallicities and the CKS. We display the CKS- Gaia 
metallicity comparison in Figure 1 and corresponding 
relations in Table 1. These homogeneous metallicities 
are available for a larger fraction of host stars, unlike 
the heterogeneous metallicities of LAMOST (Ren et al. 
2018), APOGEE (Abolfathi et al. 2018), and the CKS 
(Petigura et al. 2017) used in Berger et al. (20202). 

Table 2 shows the parameters that we input into 
isoclassify. The metallicity provenance column can 
be one of three values: Poly, Med, or None. For stars 
with Poly provenances, we used the polynomial equation 
in row 1 of Table 1 and adopt an uncertainty of 0.15 dex 
as in Berger et al. (2020a) and for those with Med prove- 
nances (outside the Tep and metallicity range of CKS 
stars), we used the linear equation in row 2 of Table 1 
and adopt an uncertainty of 0.20 dex. These uncertain- 
ties are conservative by design, as Furlan et al. (2018) 
shows that systematic metallicity uncertainties between 
spectroscopic pipelines are typically 20.1 dex. Likewise, 
the scatter in the residuals of Figure 1 is «0.11 dex. 
Stars with None provenances lack Gaia DR3 spectropho- 
tometric metallicities, and for these stars isoclassify 
defaults to a broad (20.20 dex), Gaussian solar metal- 
licity prior. 


2.3. Isochrone Fitting 


In order to determine fundamental physical parame- 
ters from the input observables adopted from Gaia DR3, 
we used isoclassify (Huber et al. 2017; Berger et al. 
2020a) in combination with a custom-interpolated PAR- 
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Figure 2. Stellar radius versus effective temperature for the 
Mann et al. (2015) M-dwarfs, where the red points represent 
the empirical data and the blue points are the isoclassify 
output using the plain PARSEC models (Bressan et al. 2012; 
Chen et al. 2019). We draw black lines connecting estimates 
for the same star. Because of this discrepancy, we use em- 
pirical corrections to constrain the properties of low mass 
stars. 


SEC evolutionary model grid (Bressan et al. 2012). We 
used kiauhoku (Claytor et al. 2020) to finely interpo- 
late the coarse models to 0.05 dex steps in metallicity 
between —2.0 to 0.5 dex, 117 masses between 0.1 and 
5.0 Mo, and ages between 0.01 Gyr and 30 Gyr. We 
used 100 equivalent evolutionary points (EEPs, Dotter 
2016) for the pre-main sequence, 250 for the main se- 
quence, 100 for the subgiant branch, and 200 for the 
red giant branch up to the tip. We used YBC bolomet- 
ric corrections (Chen et al. 2019) to compute synthetic 
photometry for our interpolated models. 

Upon applying this grid to the sample of Mann et al. 
(2015, 2019) M-dwarfs crossmatched with Gaia photom- 
etry and parallaxes, we found large discrepancies be- 
tween the measured physical parameters and those pro- 
duced by our isochrone fitting analysis. Figure 2 shows 
that large systematics are present, and that the plain 
PARSEC grid produces systematically larger and cooler 
parameters for M-dwarfs than empirically measured. 
Therefore, we used the Mann et al. (2019) M_-M_K- code, 
which uses 2MASS K, magnitudes and Gaia parallaxes 
to compute stellar masses, to compute the empirical 
masses of the Mann et al. (2015) sample. Then, we used 
those masses and the spectroscopic metallicities tabu- 
lated in Mann et al. (2015) to compute polynomial em- 
pirical relations for Gaia G, Gpp, Grp, stellar radius, 
and stellar effective temperature as a function of mass 
and metallicity. We used the Bayesian Information Cri- 
terion to determine the point at which adding additional 
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Figure 3. Empirical measurements (absolute Ggp, absolute Grp, Teg, radius) versus stellar mass colored by metallicity for 
the Mann et al. (2015) sample of M-dwarfs. We derived the masses using the Mann et al. (2019) M -M. K- code. Yellow, green, 
and indigo dashed lines show the best-fit relations for [Fe/H] = —0.6, 0.0, and 0.5 dex, and residuals show the observed values 
subtracted by the predicted value of the polynomial with a corresponding mass and metallicity, where the black dashed line 
indicates agreement. We include the numerical relations in Table 3. 
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Table 2. Gaia Kepler- TESS-Host Stellar Input Parameters 


Star ID DR3 Source ID Gsp [mag] 
kic10858832 2129500383713365248 4.5940 + 0.0030 
kic2571238 2051106987063242880 2.2374 + 0.0028 
kic8628665 2126430409811576704 5.2767 + 0.0031 
kic10328393 2130683904901068800 4.4843 + 0.0032 
epic210577548 44062348664673152 3.6642 + 0.0029 
epic220294712  2538769088655260416 2.5499 + 0.0028 
epic211711685 611385506505666688 2.7334 + 0.0028 
epic251584580  3658535369882355712 4.8086 + 0.0033 
tic429501231 3423681228084838912 3.6640 + 0.0032 
tic255685030 5501542395359045504 1.4133 + 0.0029 
tic238624131 378265951674712448 2.4537 + 0.0028 
tic391903064 5212899427468919296 9.5947 + 0.0075 


Grp [mag] c [mas] [Fe/H] [Fe/H] Prov RUWE 
13.7443 + 0.0039 1.1711 + 0.0135 0.200 + 0.150 Poly 1.003 
11.3370 + 0.0038 4.5680 + 0.0087 0.137 + 0.150 Poly 0.824 
14.4036 + 0.0039 0.9280 + 0.0185 0.072 + 0.150 Poly 1.016 
13.3366 + 0.0039 2.8866 + 0.0133 -0.003 + 0.150 Poly 1.010 
12.4527 + 0.0038 3.5165 + 0.0151 -0.027 + 0.150 Poly 0.939 
11.8115 + 0.0038 2.3343 + 0.0151 -0.092 + 0.150 Poly 1.001 
11.8419 + 0.0038 3.4644 + 0.0134 0.188 + 0.150 Poly 0.944 
13.6446 + 0.0041 2.6322 + 0.0220 -0.289 + 0.150 Poly 1.053 
12.7718 + 0.0048 1.8895 + 0.0285 -0.204 + 0.150 Poly 1.533 
10.5290 + 0.0038 6.3511 + 0.0124 0.226 + 0.150 Poly 0.990 
11.5585 + 0.0038 3.0227 + 0.0132 0.188 + 0.150 Poly 1.043 
8.7303 + 0.0061 12.6898 + 0.0105 0.154 + 0.150 Poly 0.893 


NoTE—Star ID (kic for Kepler, epic for K2, and tic for TESS), Gaia DR3 source ID, Ggp-mag, Grp-mag, parallax, metallicity, and metallicity 
provenance parameters and their uncertainties, which were provided as input to isoclassify. We also provide Gaia DR3 RUWE values for 
posterity. The [Fe/H] Prov column has three possible values: Poly, Med, and None. Poly values denote our use of the CKS-based polynomial 
method to fix the Gaia spectrophotometric metallicity (row 1 of Table 1), Med values denote our use of the median offset line (row 2 of Table 1), 
and None values denote stars that lack Gaia spectrophotometric metallicities, where isoclassify defaults to a broad, Gaussian solar metallicity 
prior. A subset of our input parameters is provided here to illustrate the form and format. The full table, in machine-readable format, can be 


found online. 


Table 3. Polynomial Empirical Relations for M-dwarfs 


Parameter Polynomial Empirical Relation 
G 24.5898 — 183.9893 M, + 1120.1562 M? — 3770.5221 M? + 7016.0366 M? — 6804.3093 M? + 2672.7206 Mê + 0.8813 [Fe/H] 
Gap 34.4834 — 318.4227 M, + 2022.7769 M? — 6932.4997 M? + 13040.8749 M1 — 12736.0413 M? + 5032.4396 M + 1.4652 [Fe/H] 
Grp 22.1215 — 165.2395 M, + 996.5293 M? — 3330.1792 M? + 6149.4185 M? — 5913.3984 M? + 2302.3600 Mê + 0.6634 [Fe/H] 
Tett 1352.4644 + 22894.5349 M, — 113739.3956 M? + 286679.2022 M? — 350226.7290 M? + 170350.2318 M? — 292.0425 [Fe/H] + 260.4901 [Fe/H]? 
Ry (0.0361 + 0.9316 M,) x (1 + 0.0638 [Fe/H]) 


NoTE—Empirical relations for M-dwarfs based on stellar mass and metallicity. The PARSEC models are defined based on their masses and metallicities, so we 
compute empirical relations for the observables as a function of mass and metallicity. The Gaia photometric relations in this table are for absolute magnitudes, 
and the Te¢e and radius relations were used to amend the log g and mean stellar densities of the models. These relations are only used for main sequence models 
with masses < 0.63 Mọ and metallicities > —0.5 dex. We plot these relations in Figure 3. We also include Gaia G for posterity. 
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Figure 4. Stellar radius versus effective temperature for 
the empirically modified, interpolated PARSEC model grid 
utilized in this work. The logarithmic number density of 
models is represented by the color of each bin. 
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Figure 5. We place a floor on Gpp-Grp color uncertain- 
ties through our empirically modified PARSEC model grid. 
Top: Teg versus Gpp-Gn» color for our model grid, with 
the logarithmic density of points illustrated by the two- 
dimensional greyscale histogram with the corresponding col- 
orbar. We plot our best-fit SmoothlyBrokenPowerLaw1D in 
red. Bottom: The red line represents the required Gpp- 
Grp minimum uncertainty to reach a 3% Te error for all 
stars, dependent on their Ggp—Grp color. 


parameters/changing the functional form no longer pro- 
duced any significant benefit. 

Figure 3 shows the Mann et al. (2015, 2019) data as a 
function of mass and colored by metallicity. In each of 
the diagrams, we do not see any strong trends in either 
the location or color of the residuals, showing that the 
mass-metallicity polynomials effectively describe the ob- 
servational data. Table 3 contains the adopted numeri- 
cal relations. For models within the mass and metallicity 
range of the Mann et al. (2015, 2019) data, we use each 
model’s mass and metallicity to replace the PARSEC 
Gap, Grp, Ter, R,, logg, and mean stellar densities 
with values derived from the empirical polynomial fits 
displayed in Figure 3 and quantified in Table 3. 

To ensure a smooth transition between the 
polynomial-modified models and the adjacent PARSEC 
models, we used a small range of masses and metallici- 
ties at which we weight the contribution of the empirical 
relations for those models. We defined this mass range 
as 0.63 < M, < 0.70 Mo, and the metallicity range as — 
0.65 < [Fe/H] < -0.50 dex, and computed the difference 
in the polynomial and evolutionary model predictions. 
We then scaled this difference based on the product of 
the mass and metallicity distance from 0.70 Mọ and 
—0.65 dex, respectively, and added this difference to 
the PARSEC model for each parameter of interest. We 
chose these particular ranges because the Mann et al. 
(2015) sample includes stars with metallicities as low 
as &:-0.6 dex and masses as large as 40.65 Mo; we 
also wanted to avoid both (1) a sharp cutoff separating 
our empirical corrections and the pre-existing models 
and (2) an extended interpolation over a larger range of 
masses/metallicities. 

Altogether, our modified PARSEC grid includes 2.8 
million models, a portion of which is represented in Fig- 
ure 4 as an HR diagram. This HR diagram includes a 
number of notable features: (1) pre-main sequence mod- 
els for low-mass stars, which can be seen as the darker 
colored dots to the right of the main sequence, (2) model 
aliasing appearing as horizontal "strips" of low-mass 
stars on the main sequence due to computational limita- 
tions of the model grid size, (3) models with empirically 
modified parameters, which appear separated from the 
evolutionary models for low-mass main sequence stars, 
and (4) the transition region from empirical relations to 
evolutionary models. 

We used isoclassify (Huber et al. 2017; Berger 
et al. 2020a) with the empirically modified PARSEC 
grid, Gaia DR3 Gep, Grp, parallaxes, positions, and 
metallicities, and the Combinedi9 map (Drimmel et al. 
2003; Marshall et al. 2006; Green et al. 2019) within 
mwdust (Bovy et al. 2016) to correct for dust extinction 
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and derive fundamental stellar parameters. Similar to 
Berger et al. (2020a), we also modified isoclassify to 
utilize uncertainties in Gpgp-Ggp that correspond to a 
~3% uncertainty in Teg by fitting a smoothly broken 
power law (astropy's SmoothlyBrokenPowerLaw1D) to 
the model grid's T;g-Gpp-Gngp relation. We chose the 
smoothly broken power law because it adequately fits 
the shape of the T;g-Gpp-Ggp curve over the range of 
Tor in the model grid with the fewest possible param- 
eters, as compared to the 12th order polynomial fit in 
Berger et al. (2020a). We then divided the Tg values by 
the derivative of this relation and multiplied by 396 to 
produce a Ggp-Grp uncertainty corresponding to 3% 
Teg uncertainty. Finally, we computed the maximum of 
the derivative Gpp-Grp uncertainty and the measured 
Gpp and Grp uncertainties added in quadrature as our 
adopted color uncertainty. Figure 5 shows this compari- 
son and the corresponding power law-derived Teg uncer- 
tainty. This ensures that we do not drastically underesti- 
mate T,g uncertainties relative to the 2.040.5% Tog sys- 
tematic uncertainty floor of interferometry (Tayar et al. 
2022) given the mmag precision of Gaia photometry and 
the imperfect power law fit to the model grid. For our 
stellar sample, we find that this procedure gives us a 
2.4% Tor uncertainty for solar-type stars. 


3. VALIDATING THE OUTPUT STELLAR 
PARAMETERS 


3.1. Accuracy of Derived Effective Temperatures, 
Radii, and Masses 


'To ensure our grid-computed stellar effective temper- 
atures and radii are accurate, we compared them to in- 
terferometric Tyg and angular diameter measurements 
for a sample of 142 stars from Huber et al. (2012), Boy- 
ajian et al. (2013), White et al. (2018), Stokholm et al. 
(2019), Rains et al. (2020), and Karovicova et al. (2020, 
2022a,b). We followed a similar procedure as above to 
crossmatch to Gaia DR3 sources and correct the pho- 
tometry and parallaxes for saturation and zeropoint dif- 
ferences, respectively. When we ran isoclassify, we 
did not correct for extinction given the interferometric 
sample's proximity to Earth. Of the 142 stars, 118 are 
within 30 pc, 131 are within 60 pc, and only one (HD 
6833) is beyond 200 pc. 

Figure 6 shows comparisons of the interferometric stel- 
lar parameters and the isoclassify-derived stellar pa- 
rameters. In general, we show agreement between the 
models and the interferometric estimates across the en- 
tire range of T,g and radius, with residual offsets of 
0.296 and 1.096 and residual scatters of 2.896 and 5.696, 
compared to the median combined interferometric and 
isoclassify uncertainties of 2.6% and 4.696, respec- 


tively. In addition, we find a systematic underestimation 
of Tag for stars hotter than 9000 K. Since the distance to 
these stars is less than 30 pc, we can rule out extinction 
as a cause for this discrepancy. Upon further inspec- 
tion, the Gaia DR3 photometry might be suspect for 
these objects, as the Gpgp-Gmgp colors are much redder 
than their effective temperature, as the current version? 
of the Pecaut & Mamajek (2013) table suggests. While 
the 0.296 and 1.096 differences between the scatter and 
median combined uncertainties are small, they suggest 
there are additional factors that could produce the ex- 
cess scatter: (1) inaccurate photometry/parallaxes, (2) 
issues with the PARSEC models, and/or (3) issues with 
the YBC synthetic photometry. 

In Figure 7 we compare our derived effective tempera- 
tures, radii, and masses to those of the Kepler LEGACY 
sample (Silva Aguirre et al. 2017), which includes esti- 
mates from a number of analysis pipelines for radii and 
masses. For Tig, radius, and mass, we find offsets of 
0.9%, 1.8%, and 1.2% and residual scatters of 1.0%, 
2.4%, and 5.0%, with differences between the individ- 
ual pipelines producing offsets and scatters of similar or 
smaller size. 


3.2. Accuracy of Derived Stellar Ages 
3.2.1. Cluster Ages 


To independently confirm our derived stellar ages, we 
used the NGC 6811 cluster data from Godoy-Rivera 
et al. (2021). We followed a similar procedure as above 
to crossmatch the NGC 6811 members found by Godoy- 
Rivera et al. (2021) to Gaia DR3 sources and perform 
corrections to the parallaxes and photometry, where nec- 
essary. We then removed stars not classified as probable 
members and those without Gaia DR3 solutions, leav- 
ing us with 338 stars. We used the Godoy-Rivera et al. 
(2021) NGC 6811 metallicity ([Fe/H] — 0.03 dex) with 
0.15 dex uncertainties for each member, and when we 
ran isoclassify, we used the mwdust allsky extinc- 
tion map as above. To ensure that our derived isochrone 
ages are reliable and minimize the potential for age con- 
taminating binaries (Berger et al. 2020a), we removed 
all stars with Tog < 6000 K and RUWE > 1.4 (Stas- 
sun & Torres 2021). For the remaining 107 stars, we 
compared stellar ages. 

Figure 8 shows our derived ages for individual stars 
and their ensemble age versus the Godoy-Rivera et al. 
(2021) NGC 6811 cluster age (0.95 Gyr). Our age dis- 
tribution is in agreement with the cluster's age. Similar 
to Berger et al. (2020a), there are a few stars at ages 


3 https: //www.pas.rochester.edu/~emamajek/EEM_dwarf_ 
UBVIJHK-colors_Teff.txt 
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Figure 6. Interferometric Tef and radius measurements versus isoclassify-derived results. Top: Radius versus Teg diagram 
showing empirical measurements in red and isoclassify results in blue, where black lines connect the separate estimates for 
the same star. Left: Direct comparison of interferometric T.g and isoclassify Tm. The black dashed line represents the 1:1 
line and the bottom panel shows the corresponding fractional residuals. Right: Same as Left but for stellar radius. 
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Figure 7. Effective temperatures, radii, and masses derived using our methodology versus those with frequency-modeled 
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The various colors/shapes represent the different pipelines. The translucent grey rounded rectangles represent the ranges of 
parameter estimates from different asteroseismic pipelines for each system, which includes not only seismic differences but also 
differences between model grids. The bottom panels show the ratio of the parameters derived in this work versus the single or 
multiple determinations of Silva Aguirre et al. (2017), which produces the curved residual structure apparent for each star in 
the right two plots. 


10 

La Isoclassify Ages 

~~ (0.9*81 Gyr) 

Nd —— NGC 6811 (0.95 Gyr) 
v 
T 0.6] 
a 
$ 
5 0.44 
(o 
T 
ac 

0.2 

0.0] TVET . 

0 1 2 3 4 5 
Age [Gyr] 


Figure 8. Age comparison for open cluster NGC 6811 at 
0.95 Gyr (Godoy-Rivera et al. 2021), located within the Ke- 
pler field. The solid red distribution represents Gaussian 
Kernel Density Estimate (KDE) of the ages of individual 
stars within each cluster as derived in this work, with the 
median and the lo confidence interval represented by the 
vertical red dashed line and shaded region, respectively. We 
use Scott’s rule (Scott 1992) bandwidths to produce the over- 
all distribution. Translucent vertical red lines represent the 
inferred ages for each star within the sample. The black, 
solid vertical line represents the cluster ages from the liter- 
ature in each panel. We only include stars with Tig > 6000 
K and RUWE < 1.4. 


>2 Gyr. After plotting the 107 NGC 6811 stars in an 
HR diagram, we find that the >2 Gyr stars are also 
the coolest and largest subgiants remaining. Although 
we already removed stars with RUWE > 1.4, some of 
these members might have very close unresolved com- 
panions (Stassun & Torres 2021) that are less massive 
and cooler, which can bias our results towards older ages 
(Berger et al. 2020a). 


3.2.2. Asteroseismic Ages 


We also compared our derived ages to those of Kepler 
stars that have asteroseismic ages as in Berger et al. 
(2020a). We utilized the “boutique” frequency-modeled 
ages from the Kepler legacy sample detailed in Silva 
Aguirre et al. (2017), which includes results from a num- 
ber of analysis pipelines. In Figure 9, the ages that we 
derive are in reasonable agreement with those provided 
by a variety of asteroseismic pipelines. The horizontal 
scatter of the colored points are typically larger than 
their reported errors, which indicates that systematic 
pipeline differences dominate. Any deviations from the 
1:1 dashed line are sufficiently accounted for by a combi- 
nation of the typical error bars (bottom right, top panel) 
and any systematic scatter depending on the asteroseis- 
mic pipeline one chooses. The structure in the residuals 
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Figure 9. Ages derived using our methodology versus those 
with frequency-modeled asteroseismic ages from a variety of 
pipelines (Silva Aguirre et al. 2017). The black dashed line 
represents agreement. The various colors/shapes represent 
the different pipelines. The translucent grey rounded rect- 
angles represent the ranges of age estimates from different 
asteroseismic pipelines for each system, which includes not 
only seismic differences but also differences between model 
grids. The bottom panel is the ratio of the age derived in 
this work versus the multiple determinations of Silva Aguirre 
et al. (2017), which produces the curved residual structure 
apparent for each star. We have also plotted median error 
bars in the right-hand portion of the top panel, where, from 
bottom to top, the error bars represent the median uncer- 
tainties of stars with isochrone ages between 0-4 Gyr, 4-8 
Gyr, and >8 Gyr, respectively. 


in the bottom panel results from their representation as 
a ratio between the age derived in this work and those 
of Silva Aguirre et al. (2017). In addition, the astero- 
seismic ages do not fall above the age of the Universe 
(likely due to a model grid age-cutoff). We also see that 
at intermediate ages (4-8 Gyr), we may slightly under- 
estimate stellar ages, but our solutions are well-within 
our reported uncertainties. 

Ultimately, we report a median offset and scatter of 
0.796 and 3196 between our isochrone-derived ages and 
the asteroseismic ages of Silva Aguirre et al. (2017), re- 
spectively. We conclude that our isochrone-derived ages 
are consistent with ages determined through more pre- 
cise methods within the uncertainties that we report. 
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4. RESULTS 
4.1. Stellar Properties 


We present our homogeneous stellar properties in Ta- 
ble 4 and Figure 10. Figure 10 shows the host star pop- 
ulations broken up into separate panels for each mission 
and the combined host star sample in the bottom right. 
Most Kepler host stars are solar-like stars, as was the 
priority of the Kepler mission's target selection (Borucki 
et al. 2010; Batalha et al. 2010) and as has been con- 
firmed more recently thanks to Gaia photometry and 
parallaxes (Wolniewicz et al. 2021). Of the 3248 Ke- 
pler host stars, only 17 are hotter than 8000 K, while 95 
are cooler than 4000 K. Relative to the K2 and TESS 
host stars, there are more giants with planet candidates 
in the Kepler host star sample. This is likely because 
of Kepler’s superior sensitivity and baseline when com- 
pared to the shorter continuous baselines and increased 
systematic noise/smaller apertures of K2/ TESS. 

The upper right panel of Figure 10 displays the K2 
host star sample. There are fewer K2 hosts than both 
Kepler and TESS, mostly because we used the homo- 
geneously derived input sample from Zink et al. (2021). 
Due to differences in target selection, most K2 host stars 
are lower mass than those observed by Kepler. Eighty- 
nine K2 stars are cooler than 4000 K and 216 are cooler 
than 5000 K, as K2 target selection was community- 
led and prioritized K and M-dwarfs in its search for 
planets (Dressing & Charbonneau 2013; Howell et al. 
2014; Dressing & Charbonneau 2015; Huber et al. 2016b; 
Cloutier & Menou 2020). Only one star in our K2 sam- 
ple is hotter than 8000 K. In the lower left panel, we 
plot the T'ESS host star distribution. T'ESS host stars, 
perhaps unsurprisingly, contain the largest diversity of 
stellar types, given T'ESS's all-sky observations. Host 
stars range in Tog from 2934 K up to 18740 K, and of 
the 4180 hosts, 317 are cooler than 4000 K and 179 are 
hotter than 8000 K. However, we caution that a large 
fraction (78596) of the TESS host stars are hosts to 
planet candidates and hence have not been subject to 
the scrutiny of ground-based follow-up or comprehen- 
sive vetting routines. 

Together, our Kepler4-K24- TESS host star sample, 
displayed in the bottom right panel of Figure 10, repre- 
sents both the largest and most diverse transiting exo- 
planet host star sample with homogeneous fundamental 
parameters yet. This enables a direct comparison of 
the host stars and their stellar properties and the con- 
struction of larger host samples across a wide range of 
parameter space. For instance, we count 492 host stars 
with Tee > 6000 K and isochrone ages «1 Gyr. These 
stars represent an interesting sample for further study, 


although we note the potential for many of the planet 
candidates orbiting these stars, especially those > 10000 
K, to be false positives. The ages of exoplanet host stars 
are particularly useful for determining how exoplanets 
evolve over time (Soderblom 2010; Mann et al. 2017; 
Berger et al. 2018a, 2020b, 2022; David et al. 2021). 

We plot the distribution of masses and ages in Fig- 
ure 11. As is clear in Figure 10, the mass (and age) 
ranges vary from sample to sample, where Kepler host 
stars are mostly close to one solar mass, while K2 hosts 
shift to lower masses and TESS hosts include a broader 
distribution of stellar masses. In stellar age, we see that 
the Kepler and K2 distributions are similar except for 
the peak in stellar age between 13 and 15 Gyr - this is 
because K2 focused more heavily on low-mass K and M- 
dwarfs, which do not produce meaningful isochrone ages 
and tend to have ages roughly half the maximum age of 
our model grid (30 Gyr). The TESS histogram displays 
a large number of hosts younger than 3 Gyr, most of 
which are stars in the high mass tail of the TESS mass 
histogram. 


4.2. Exoplanet Properties 


As in Berger et al. (2020b), we utilized the Cumulative 
KOI table* at the NASA Exoplanet Archive (Batalha 
et al. 2013; Burke et al. 2014; Rowe et al. 2015; Mul- 
lally et al. 2015; Coughlin et al. 2016; Thompson et al. 
2018; NASA Exoplanet Archive 2021) catalog-provided 
R,/R, and orbital period values, where available, in 
combination with our stellar parameters and their uncer- 
tainties to compute planet radii, semimajor axes, and in- 
cident fluxes. Both the Cumulative KOI and Zink et al. 
(2021) catalogs provide R,/R, values for the Kepler and 
K2 planet samples, but the T'ESS Project Candidate ta- 
ble (e.g. Guerrero et al. 2021) at the NASA Exoplanet 
Archive does not. Therefore, we used this TESS tables? 
stellar and planet radii to compute R,/R, which we then 
multiplied by our stellar radii to compute TESS planet 
radii We do not use transit depth measurements to 
compute T'ESS planet radii because transit depth mea- 
surements do not directly translate into planetary radii 
for planets that have large impact parameters. 

We present our planet properties in Table 5 and Fig- 
ure 12. Similar to the stellar populations in 84.1, the 
planet populations in the planet radius-incident flux dia- 
grams differ from panel-to-panel. In the upper left panel 
containing Kepler planets, we can see a few, clear sub- 
populations of planets. We see two populations of plan- 
ets smaller than 5 Ra, separated at by the planet radius 


^ accessed 11/5/21 
5 accessed 7/27/22 
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Figure 10. Stellar radius versus effective temperature for the Kepler, K2, TESS, and combined host star samples. 


gap at ~2 Ra: these are the super-Earths (<2 Re) and 
sub-Neptunes (22 Rgand <4 Re), as identified by Ful- 
ton et al. (2017). There are many super-Earths due to 
Kepler’s sensitivity and baseline, and the sub-Neptune 
population appears to be well-defined given the drop- 
off in planets at >5 Rg. Finally, we observe a group 
of Jupiter-sized planets (10-20 Rg) at a wide range of 
incident fluxes. 

We show the K2 planet population in the upper right 
panel of Figure 12. There are far fewer K2 planets 
than in either Kepler or TESS, and while we see the 
sub-Neptune population, K2 did not detect many super- 
Earths. Hence, the planet radius gap is not clear in the 
K2 sample due to the small number of super-Earths, al- 
though it was detected in Hardegree-Ullman et al. (2020) 


and Zink et al. (2021). The differences in planet radii 
must come from the utilized stellar radii; we used the 
same R,/R, values as Zink et al. (2021) but used a 
predominantly photometric isochrone fitting approach 
as compared to the spectroscopically-dependent random 
forest classifier of Hardegree-Ullman et al. (2020). Com- 
paring the two sets of stellar radii, we find a median off- 
set of x3% and typical scatter of x5% with some larger 
outliers. The observation that the K2 planet radius gap 
can disappear because of a statistically consistent yet 
slightly different set of stellar parameters suggests it is 
not yet robustly determined, unlike the Kepler planet 
radius gap. 

We plot the T'ESS planets in the bottom left panel, 
and while there are more TESS planets than Ke- 
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Figure 11. Histograms of masses and ages for the Kepler, K2, and TESS host star samples, where the left column contains 
stellar masses and the right column contains stellar ages. The Kepler, K2, and TESS, samples are displayed in each row and 


ordered, respectively, from top to bottom. 


pler planets, Kepler has more super-Earths and sub- 
Neptunes. The vast majority of TESS planets are 
hot Jupiters, which appear to become more inflated at 
higher incident flux (Miller & Fortney 2011; Thorngren 
et al. 2016; Grunblatt et al. 2016, 2017, 2019), although 
TESS has also detected quite a few sub-Neptunes. Due 
to its shorter observation baselines and reduced sensitiv- 
ity relative to Kepler (Ricker et al. 2014), TESS has not 
detected as many super-Earths. Hence, we do not see 
a clear gap separating super-Earths and sub-Neptunes. 
One explanation for the lack of a clear gap in both the 
K2 and TESS samples are their less-precise Ry / R, val- 
ues (Petigura 2020). As above, we caution that a large 
fraction of the T'ESS planets are planet candidates and 
upon follow-up some number will be identified as false- 
positives. 

'The combined planet sample in the bottom right panel 
of Figure 12 represents the largest transiting exoplanet 


sample with homogeneously derived planet parameters 
yet. This enables us to increase the number of planets 
in notable regions of parameter space, such as the hot 
sub-Neptunian desert (Rp = 2.2-3.8 Re and Fp > 650 
Fg, Lundkvist et al. 2016; Berger et al. 2018b, 2020b), 
the habitable zone (Fp = 0.25-1.50 Fg, Kane et al. 2016; 
Berger et al. 2018b, 2020b), and the planet radius gap 
(Fulton et al. 2017; Fulton & Petigura 2018; Van Eylen 
et al. 2018; Berger et al. 2018b, 2020b; Petigura 2020; 
Petigura et al. 2022). Moreover, the populations of in- 
flated Jupiters, both hot and warm/cold (Berger et al. 
2018b, 2020b), warrant additional scrutiny. 

We take a closer look at the Kepler, K2, TESS, and 
combined planet radius histograms in Figure 13. Sim- 
ilar to our conclusions based on the various panels in 
Figure 12, we see no clear gap in either the K2 or TESS 
planet radius histograms. We emphasize that these are 
raw numbers and are not completeness-corrected. Com- 
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Figure 12. Planet radius versus incident flux for the Kepler, K2, and TESS, and combined planet samples. While most planets 
are shown, there are 1 Kepler and 68 TESS planets that experience incident fluxes < 0.01 Fẹ. 


paring the Kepler and combined samples, it is clear that 
Kepler provides superior gap resolution, and the addi- 
tion of K2 and TESS planets makes it less-clear. 


5. CONCLUSION 


We presented the Gaia-Kepler-TESS-Host Stellar 
Properties Catalog and the corresponding catalog 
of homogeneous exoplanet properties, which contain 
7993 stars and 9324 planets, respectively. We used 
isoclassify and Gaia DR3 to compute precise, ho- 
mogeneous Teg, logg, masses, radii, mean stellar den- 
sities, luminosities, ages, distances, and V-band extinc- 
tions for 3248, 565, and 4180 Kepler, K2, and TESS 
stars, respectively. We compared our stellar properties 


to studies using fundamental and precise constraints, 
such as interferometry and asteroseismology, and find 
good agreement between our Tyg, radii, and ages and 
those in the literature. In addition, we provide planet 
radii, semimajor axes, and incident fluxes for 4281, 676, 
and 4367 Kepler, K2, and TESS planets, respectively, 
and find that the exoplanet radius gap is less prominent 
in K2 and TESS and combined samples than it is in the 
Kepler sample alone. We provide our stellar and planet 
catalogs as machine-readable tables associated with this 
article and as a High-Level Science Product (HLSP) at 
the Mikulski Archive for Space Telescopes (MAST). Fi- 
nally, we identify over 1000 hot Jupiter planet candi- 
dates, 150 planets within the hot sub-Neptunian desert, 
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Figure 13. Planet radius histograms for the Kepler, K2, 
and TESS, and combined planet samples. 


and more than 400 young host stars as potential oppor- 
tunities for testing the theories of planet formation and 
evolution. 


Software: Astropy (Astropy Collaboration et al. 
2013, 2018, 2022), emcee (Foreman-Mackey et al. 2013), 
GNU Parallel (Tange 2022), kiauhoku (Claytor et al. 
2020), matplotlib (Hunter 2007), numpy (Harris et al. 
2020), pandas (McKinney 2010), scipy (Virtanen et al. 
2020), TOPCAT (Taylor 2005) 

Facility: Exoplanet Archive 
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Table 4. Gaia-Kepler- TESS-Host Stellar Output Parameters 


Star ID Tere [K] log g [dex] [Fe/H] M, [Mo] R, [Ro] Px [Po] L, [Lo] Age [Gyr] Distance [pc] Ay [mag] 
kic10858832 5812109 44241003 — 0.140*0 131 — 1.054*0-022  1.040*0 05)  0.927*0 109 — 1322*0047 — 2.86!2 72 853.5*19$ 0.110 
kic2571238 5444+990  4.510*0022 — 0.055*0 138 — 0.915*0 084 — O.880*0 026 — 1.334*0 41, 0.62070 019 — 4.29* 827 218.9* 21 0.000 
kic8628665 5895*115  4.450*0020  .0.021*0 139  1.034*0 029  1.000*'0 059  1.024*0 102 1.00370 002 — 2.31*7 23 — 1079.5728 2 0.275 
kic10328393 4818***.  4.5906*0 022  -0.002*0-120  0.735*0 039  0.715*0 024  2.001*022,  0.250*0005  *T.98*1l1? 346.7153 0.027 
epic210577548 5381*111 4.53170 064  -0.008*0:129 — 0.854*0 0:0 0.829*0 055 — 1.486'0 250  0.525*0 01a — 5.013 101 284.3*13 0.604 
epic220294712 62947735  4.360*0 026 — -0.1277* 1482-0955 LUGO Ebay Oris pete 1:915 EG OSs :88*7:57. 428.0* 6.6 0.192 
epic211711685  5484*95,  4.496*0023  OQ.111*0 l3.  0.936*0 040  0.906*00260 — 1.256'0 159 0.676*0 02. = 4.00* $05 288.7512 0.000 
epic251584580  4830*1353  4,.601*0044  .9,352*0.128 — 9,696*0.021 — 9,68g*0.028  2109*0304  Q333*0014  *g5g*lL13 ^ grgg*62 0.137 
ticd29501231 — 5779*102  4.371*0 05.  -0.186*0 125  0.914*0 040  1.025*0 049  0.828*0128 1.07370 095 8.014869% 527.8* 111 0.110 
tic255685030 = 5519*91  4.502*0025  0.120*'0 134  0.957*0 044. 0.9009*0 02,  1.270'0 114  0.696'0 02; — 3167254 157.5*2:2 0.053 
110238624131 — 5924*11?  4.384*0 022 — 0.148*0 138 — 1.006*0 065  1.110*'0 055 0.7920 102 1.37670 022 — 2.84*2 12 330.6* 12 0.247 
tic391903064  5721*197  4.451*0 029 — 0.088*0 135 — 1.012*0 000 0.987 *0 028 1.03870 755 0.95078 025 — 3.30* 221 78.8*11 0.089 

NoTE—KIC ID, effective temperature, surface gravity, surface metallicity, stellar mass, stellar radius, density, luminosity, age, distance, and V-magnitude 
extinction, output from our isochrone fitting routine detailed in 82. Asterisks are appended to stellar ages that are uninformative (Tere « 5000 K). A subset 
of our output parameters is provided here to illustrate the form and format. The full table, in machine-readable format, can be found online. 


Table 5. Gaia-Kepler- TESS Planet Parameters 


Star ID Planet ID Disposition P [days] Ry/ R, Ry [Re] a [AU] Fp [Fẹ] 
kic10797460 K00752.01 CP 9.488036* 0-000028 0.02234 0 a3 2.201012 0.085775 D-00119 99.190* 7: 220 
kic10797460 K00752.02 CP 54.418383 *0 000248 0.02795 70 00908 2.761090  Q.27482* 0 bees 9.661*0 74 
kic10811496 K00753.01 PC 19.899140*0 000015 — 0.15405* 5 012.&  18.827015.0? — 0.146287 0 00209 61.417* $ 270 
kic10854555 K00755.01 CP 2.525592* 0 000004 0.024067 0 barre 2.42*0-10 — 0.035507 0 00095 629.842 10-222 
epic201155177 — epic201155177.01 KP 6.689040*0-000799 — 0.02850 0 Qo100 2.151697 . 0.06105 000109 49.117* 555 
epic201208431 — epic201208431.01 KP 10.004208* 9000221 — 0.03520" 0 00100 DATE ne 0077375 Dose d7:975 0059 
epic201247497 — epic201247497.01 PC 2.754012 0 000122. — 0.091707 0 02500 6.82*?.09 — 0.03339" 0 00050 116.013 * 11 210 
epic201295312 — epic201295312.01 KP 5.656766" 0 000278 0.017101000090 2.781017 0.065230 00100 567.312* 36 121 
tic88863718 toil001.01 PC 1.931671*0 000009 — 0.04701*0 01480 0,055220 ./0:0835170 00023. “T8087 77 tae, ose 
tic65212867 toil007.01 PC 6.998921* 0000014 0.050124000351 14.40*122  0.08440*0 00097 — 1506.786* 2. 120 
tic231663901 toil01.01 KP 1.430369*0-000001 — 0.13623*0 00954 13.23*109  0.02429*0 00021 — 1301.559*55 105 
tic114018671 toil011.01 PC 2.469696* 0 000279 — 0.014157 0 bi gee 1.381785 —0.03468*0 00092 547.7407 22 957 


NoTE—Host star ID, planet ID, disposition, orbital period, Ry / R,, planet radius, semimajor axis, and incident flux values and their uncertainties 
for all planets contained within our catalog. Dispositions are designated as KP, CP, PC, APC, or LPPC for Known Planet, Confirmed Planet, 
Planetary Candidate, Ambiguous Planetary Candidate, or Long Period Planetary Candidate, respectively. A subset of our planet parameters is 
provided here to illustrate the form and format. The full table, in machine-readable format, can be found online. 
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